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Abstract. We present recent results on leptonic B decays using data collected by the BaBar detector 
at the PEP-II asymmetric-energy e+e- collider at the Stanford Linear Accelerator Center. We 
report searches for the B'^ — > r^i/ decay based on two statistically independent data samples. 
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1 Introduction 

The purely leptonic decay — > r+i/ [1] is sensitive to 
the product of the B meson decay constant Jb , and the 
absolute value of Cabibbo-Kobayashi-Maskawa matrix 
element Kb [2, 3]. In the Standard Model (SM), the de- 
cay proceeds via quark annihilation into a boson, 
with a branching fraction given by: 
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where Gp is the Fermi constant, Tg+ is the B~^ life- 
time, and rrig and rrir are the B~^ meson and r lepton 
masses. 

The process t^v is also sensitive to exten- 

sions of the SM. For instance, in two-Higgs doublet 
models [4] and in the MSSM [5, 6] it could be medi- 
ated by charged Higgs bosons. The branching fraction 
measurement can therefore also be used to constrain 
the parameter space of extensions to the SM. The Belle 
Collaboration has reported evidence from a search for 
this decay and the branching fraction was measured to 
be (1.79t|]:^^(stat.)l:[J-j^(syst.)) x lO"'^ [7]. We present 
here two results from the BABAR collaboration using 
a sample of 383 x 10*^ 7(45*) -> BB decays, based 
on the reconstruction of a semileptonic [8] and of an 
hadronic [9] B decay on the tag side. 



2 Tag B Reconstruction 

In the hadronic tags analysis, the tag B candidate 
is reconstructed in hadronic B decay modes B~ —^ 
£)(*)o^-^ where denotes a system of charged and 
neutral hadrons with total charge —1 composed of 
niTT^, n2K^, n^K^, 71477°, where ni 4-712 < 5, < 2, 
and 714 < 2. We reconstruct D*° D°Tr°,D°-f; D° 
K-Tr+, K-n+Tr°, K-n+TT-n+, K^tt+tt- and X° 



TT+vr . Tag B candidates are required to be kinemat- 
ically consistent with decay from an T{AS) using the 
beam energy-substituted mass rriEs ~ \/s/4: — Pg and 
the energy difference AE = Eb — a/s/2. Here y/s is the 
total energy in the T(4S') center-of-mass (CM) frame, 
and pb and Eb denote, respectively, the momentum 
and energy of the tag B candidate in the CM frame. 
The purity V of each reconstructed B decay mode is 
estimated, using on-resonance data, as the ratio of the 
number of peaking events with tmes > 5.27 GeV/c^ 
to the total number of events in the same range. If 
multiple tag B candidates are reconstructed, the one 
with the highest purity V is selected. If more than one 
candidate with the same purity is reconstructed, the 
one with the lowest value of \AE\ is selected. The set 
of decay modes used is defined by the requirement that 
the purity of the resulting sample is not less than 30%. 

The background consists of e'^e~ — *■ qq events and 
other T{4:S) B^B° or B+B- decays in which the 
tag B candidate is mis-reconstructed using particles 
coming from both B mesons in the event. To reduce 
the e+e" qq background, we require jcos^^^^l < 
0.9, where O'^g is the angle in the CM frame between 
the thrust axis [10] of the tag B candidate and the 
thrust axis of the remaining reconstructed charged and 
neutral candidates. 

In order to determine the number of correctly re- 
constructed B^ decays, the background events are_clas- 
sified in four categories: e+e" — > cc; e+e~ uu, dd, ss; 
r(45) ^ B°B°; and T(4S') ^ B+B^ . The tties shapes 
of these background distributions are taken from MC 
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e+e^ uu, dd, ss backgrounds is taken from off- 
resonance data, scaled by the luminosity and corrected 
for the different selection efficiencies evaluated with 
the MC. The normalization of the B°S°, B+B' com- 
ponents are obtained by means of a fit to the 771^3 
distribution in the data sideband region (5.22 GeV/c^ < 
niBS < 5.26 GeV/c^). The number of background events 
in the signal region {m-^s > 5.27GeV/c^) is extrapo- 
lated from the fit and subtracted from the data. We 
estimate the total number of tagged B's in the data 
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to be Nb = (5.92 ± O.ll(stat)) x 10^ Figure [T] shows 
the tag B candidate rngs distribution, with the combi- 
natorial background, estimated as the sum of the four 
components described above, overlaid. 
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Fig. 1. Distribution of the energy substituted mass, rriES, 
of the tag B candidates in data. The combinatorial back- 
ground is overlaid. 



In the semileptonic tag analysis, the tag B is re- 
constructed in semileptonic B decay modes B^ — s- 
D^£~P(X, where i is either an electron or a muon, 
and X can be either nothing, a tt" or a photon. Events 
where the best tag candidate is consistent with neutral 
B decay are rejected. The identified electron or muon 
in the D'^i candidates are required to have momen- 
tum above 0.8 GeV/c in the e+e~ center-of-mass (CM) 
frame. The flight direction of the D° is required to in- 
tersect with the lepton track. We reconstruct the 
candidates in four decay modes: K~Tr~^ , K~ tt^ i:^ t:^ , 
i^~7r+7r'^, and /CgTr+Tr", only considering candi- 
dates decaying to charged pions. The 7r° candidates 
are required to have invariant masses between 0.115 
and 0.150 GeV/c^ and the photon daughter candidates 
of the tt'^ must have a minimum laboratory energy of 
30MeV and have shower shapes consistent with elec- 
tromagnetic showers. The mass of the reconstructed 
13° candidates in the K~tt'^ , K~ tt'^'k~ tt'^ , and K^^tt^tt^ 
modes is required to be within 20 MeV/ of the nom- 
inal mass [11], while in the K^tt^-k^ decay mode the 
mass is required to be within 35 MeV/c^ of the nominal 
mass. Furthermore, the sum of the charges of all the 
particles in the event must be equal to zero. 

We calculate the cosine of the angle between the 
D^l candidate and the B meson as 
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constructed the D'^f candidate with the largest prob- 
ability of originating from a single vertex is selected. 

From signal MC we estimate the tag reconstruction 
efficiency to be (6.64 ± 0.03) x lO'^, where the error 
is due to the statistics of the signal MC sample. This 
corresponds to a tag B yield of (2.54 ± 0.03) x 10^. 



3 Selection of Signal Events 

After the reconstruction of the tag B meson, the rest 
of the event (recoil) is examined for 5+ t^v decays. 
We require the presence of only one well-reconstructed 
charged track (signal track) with charge opposite to 
that of the tag B. The signal track is required to have 
at least 12 hits in the drift chamber, momentum trans- 
verse to the beam axis, pr, greater than O.lGeV/c, 
and the point of closest approach to the interaction 
point less than 10 cm along the beam axis and less 
than 1.5 cm transverse to it. 

The T lepton is identified in four decay modes con- 
stituting approximately 71% of the total r decay width: 



T"*" e^vv, r"*" — > pL^vv, t"*" — > ■n'^v^ and r"*" 



signal 



TT+TT F. Particle identification criteria on the 
track are used to separate the four categories. 

The T+ — > 7r+7r°I7 sample is obtained by associ- 
ating the signal track, identified as pion, with a 7r° 
reconstructed from a pair of neutral clusters with in- 
variant mass between 0.115 and 0.155 GeV/c^. In the 
hadronic tag analysis and tt+tt*^ energy is required to 
be greater than 250 MeV and in case of multiple candi- 
dates, the one with largest center-of-mass momentum 
p*+.^o is chosen. In the semileptonic tag analysis the 
energy of the neutral clusters is required to be greater 
than 50 MeV. 

The selection is further refined with additional re- 
quirements exploiting the kinematics of the signal, in- 
cluding cuts on the momentum of visible r"*" decay 
products, the missing momentum, the missing mass, 
the invariant mass of the tt+tt" candidate, charged 
tracks and tt" candidates multiplicities. The semilep- 
tonic tag analysis also includes a veto on signal of 
in the calorimeter and in the muon detector. 

We optimize the selection by maximizing s/V s + h 
using the MC simulation, where h is the expected back- 
ground and s is the expected number of signal events 
in the hypothesis of a branching fraction of 1 x lO^''. 
The optimization is performed separately for each r 
decay mode and with all the cuts applied simultane- 
ously in order to take into account any correlations 
among the discriminating variables. The semileptonic 
tag analysis uses the PRIM algorithm [12] to find the 
optimal set of cuts. 

For both the analyses, the most powerful discrimi- 
nating variable is -Eoxtra defined as the sum of the ener- 
gies of the neutral clusters not associated with the tag 
B or with the signal 7r° from the r"*" tt+tt^F mode. 



where {E]jo£, Pd^i) is the four-momentum of the D''^i 
candidate in the CM frame, and rrijjoi and ms are the 
invariant masse of the D'^i candidate and the B'^ me- 
son nominal mass [11], respectively. We expect cos 6'B_£)o^and passing a minimum energy requirement. In the 
for correctly reconstructed tag B candidates to be in hadronic tag analysis the required energy depends on 
the range [—1, 1], whereas combinatorial backgrounds the selected signal mode and on the calorimeter region 
can have values outside this range. We select events involved and varies from 50 to 70 MeV. In the semilep- 
with —2.0 < cos9 B_j)o^ < 1.1, If multiple tag are re- tonic tag analysis, the minimum energy requirement 
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is fixed to 20 MeV, and the energy of charged tracks 
not associated with the tag B or the signal candidate 
are included in the -Boxtra computation. Signal events 
peak at low -Eoxtra values, whereas background events, 
which contain additional sources of neutral clusters, 
are distributed toward higher E'cxtra values. 

The total selection efficiency is estimated from sig- 
nal MC to be e = (9.8 ± 0.3)%, and e = (12.7 ± 0.2)% 
for the hadronic tag analysis and the semileptonic tag 
analysis, respectively. 



Table 1. Observed number of on-resonance data events 
in the signal region compared with the number of expected 
background events, for the hadronic tag analysis. 



r decay mode 


Expected background 


Observed 


r"*" — * e^vV 


1.47 ± 1.37 


4 




1.78 ± 0.97 


5 


r"'" — > tt'^V 


6.79 ± 2.11 


10 




4.23 ± 1.39 


5 


All modes 


14.27 ± 3.03 


24 



4 Background yield 



In the hadronic tag analysis, we define a sideband 
region 0.4 GeV < -Eextra < 2.4 GeV and a r mode- 
dependent signal regions. We perform an extended un- 
binned maximum likelihood fit to the toes distribution 
in the i?cxtra data sideband region of the final sam- 
ple. For the peaking component of the background we 
use a probability density function (PDF) which is a 
Gaussian function joined to an exponential tail (Crys- 
tal Ball function) [13]. As a PDF for the non-peaking 
component, we use a phase space motivated threshold 
function (ARGUS function) [14]. From this fit, we de- 
termine a peaking yield A'^j^, ^ ^ and signal shape pa- 
rameters, to be used in later fits. The same procedure 
is applied to B^B~ MC events which pass the final 
selection to determine the peaking yield N^^°"^'^ . To 
determine the MC peaking yield in the -Ecxtra signal 
region iV^Jf'^^^, we fit toes in the Ec^tra. signal region 
of the B~^B~ MC sample with the Crystal Ball pa- 
rameters fixed to the values determined in the -Ecxtra 
sideband fits described above. Analogously, we fit the 
Toes distribution of data in the -Ecxtra signal region 
to extract the combinatorial background ricomb, evalu- 
ated as the integral of the ARGUS shaped component 
in the toes > 5.27 GeV/ c'^ region. The total expected 
background in the signal region is determined as 
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In the semileptonic tag analysis, the sideband region is 
defined by i?cxtra > 0.5 GeV, and the signal regions in 
B-oxtra are signal mode-dependent. Using the number of 
events in the sideband (A^MC.sb) and signal (A^MC.sig) 
regions from MC simulation and the number of data 
events in the sidebands A^data.sb , we estimate the num- 
ber of expected background events in the signal region 
in data A^n^n siir 



exp,sig 



MC,sig 
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The background estimate is validated using side- 
bands in the D'^ mass distribution. 



5 Results 



We measure the yield of events in each decay mode in 
on-resonance data. Table [1] reports, for the hadronic 
tag analysis, the number of observed events together 
with the expected number of background events, for 
each r decay mode. Figure [2] shows the i?cxtra distri- 
bution for data and expected background at the end of 
the selection. The signal MC, normalized to a branch- 
ing fraction of 3 x 10~^ for illustrative purposes, is 
overlaid for comparison. The i?oxtra distribution is also 
plotted separately for each t decay mode. 
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Fig. 2. -Bcxtra distribution after all selection criteria have 
been applied. The on-resonance data (black dots) distri- 
bution is compared with the total background prediction 
(continuous histogram). The hatched histrogram repre- 
sents the combinatorial background component. —* 
t'^v signal MC (dashed histogram), normalized to a 
branching fraction of 3 x 10"'' for illustrative purposes, 
is shown for comparison. 



Table [2] lists, for the semileptonic tag analysis, the 
number of observed events in on-resonance data in the 
signal region, together with the expected number of 
background events in the signal region. Figure [3] shows 
the i?oxtra distribution for all data and MC in the sig- 
nal region, with signal MC shown for comparison. 

In both the analyses, we combine the results on 
the observed number of events ni and on the expected 
background hi from each of the four signal decay modes 
{rich) using the estimator Q — L{s + b) / L{b), where 
C{s + b) and C{b) are the likelihood functions for sig- 
nal plus background and background-only hypotheses. 
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Fig. 3. -Bcxtra distribution. All selectiorffti'teria liave been 
applied and all signal modes combined. Background MC 
(solid histogram) has been normalized to the luminosity of 
the on-resonance data (black dots), and then additionally 
scaled according to the ratio of predicted background from 
data and MC. B'^ t^v signal MC (dotted histogram) is 
normalized to a branching fraction of 10"'^ and shown for 
comparison. 

Table 2. Observed number of on-resonance data events 
in the signal region are shown, together with number of 
expected background events, for the semileptonic tag anal- 
ysis. 



r 


Expected background 


Observed 


r"*" — > e^vV 


44.3 ± 5.2 


59 


r"*" — > fi'^vV 


39.8 ± 4.4 


43 


t'^ — > tt'^V 


120.3 ± 10.2 


125 


r+ ^ Tv+7v°V 


17.3 ± 3.3 


18 


All modes 


221.7 ± 12.7 


245 



In the hadronic tag analysis, we measure the branch- 
ing fraction 

B{B+ ^ T+v) = {l.Stol ± 0.4 ± 0.2) X 10-^ (7) 

where the first error is statistical, the second is due 
to the background uncertainty, and the third is due 
to other systematic sources. Taking into account the 
uncertainty on the expected background, as described 
above, we obtain a significance of 2.2 cr. 

In the semileptonic tag analysis, we measure the 
branching fraction 

B{B+ T+v) = (0.9 ± 0.6 ± 0.1) X 10"^ (8) 

where the first error is statistical, and the second is 
due to the systematics uncertainties. 

The combination of the two results yields: 

B{B+ -> T+v) = (1.2 ± 0.4 ± 0.3 ± 0.2) x 10"^ (9) 

where the first uncertainty is statistical, the second 
is the uncertanty on the expected background and the 
third is from the other source of systematic effects. The 
significance of the combined result is 2.6 a including 
the uncertainty on the expected background (3.2 a if 
this uncertainty is not included). 
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respectively: 
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(5) 

The estimated number of signal candidates Si in data, 
for each decay mode, is related to the JS"*" t^v 
branching fraction by: 
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where N is the number of tag mesons correctly 
reconstructed, e"^^ and £*fg are the tag B efficiencies 
in generic BB and signal events respectively, and Si 
are the signal selection efficiencies for each channel, 
including the t+ branching fractions. In the hadronic 
tag analysis, we find the ratio from MC simulation to 
be elj/es^ = 0.939 ± 0.007(stat.). 

We estimate the branching fraction (including sta- 
tistical uncertainty and uncertainty from the back- 
ground [15]) by scanning over signal branching frac- 
tion hypotheses and computing the value of L{s -\- 
h)/C{h) for each hypothesis. The branching fraction 
is the hypothesis which minimizes the likelihood ra- 
tio — 21nQ = — 21n(£(s -I- b)/C{b)), and we determine 
the statistical uncertainty by finding the points on the 
likelihood scan that occur at one unit above the min- 
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